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ABSTRACT 

The unusual properties of the bursting X-ray pulsar GRO J1744-28 are 
explained in terms of a low-mass X-ray binary system consisting of an evolved 
stellar companion transferring mass through Roche-lobe overflow onto a neutron 
star, implying that the inclination of the system is <^ 18°. Interpretation of 
the QPO at frequency z/qpo = 40 Hz using the beat-frequency model of Alpar 
& Shaham and the measured period derivative P with the Ghosh & Lamb 
accretion-torque model implies that the persistent X-ray luminosity of the 
source is approximately equal to the Eddington luminosity and that the neutron 
star has a surface equatorial magnetic field B eq = 2 x 10 10 (z^qpo/40 Hz) -1 G 
for standard neutron star parameters. This implies a distance to GRO J1744-28 
of = 5 (Vqpo/40 Hz) 1 / 6 ^ 1 / 2 kpc, where b < 1 is a correction factor that depends 
on the orientation of the neutron star. 



Subject headings: accretion — stars:binary:close - - stars:neutron 
pulsars:individual (GRO J1744-28) — X-rays:bursts 



Introduction 



GRO J1744-28, a 2.1 Hz X-ray pulsar with unusual bursting properties, was discovered 
on 2 December 1995 in the direction of the Galactic Center with the Burst and Transient 
Source Experiment (BATSE) on the Compton Gamma Ray Observatory (Kouveliotou et 
al. 1996a; Finger et al. 1996a). The source initially bursted about every 3 minutes, then 
declined in frequency to about 30 per day for the first month, after which it increased in 
frequency to about 40 per day in the period 8-18 January 1996 (Fishman et al. 1996). The 
Rossi X-ray Timing Explorer (RXTE) (Swank et al. 1996) localized the source to galactic 
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coordinates I = +0.02°, b = +0.3°. Both the persistent and bursting X-ray emissions show 
pulsations with a 0.467 s period (Kouveliotou et al. 1996b) and have similar hard X-ray and 
soft gamma-ray spectra (Briggs et al. 1996; Strickman et al. 1996) with a color temperature 
~ 10 keV, characteristic of X-ray pulsars. 

In this Letter, we show that the observations are best explained by a neutron star 
accreting matter through Roche lobe overflow from an evolved low-mass stellar companion. 
We derive values for the magnetic field stength and accretion luminosity by noting that the 
beat frequency model (Alpar & Shaham 1985) for QPOs and the accretion-torque model 
(Ghosh & Lamb 1979) for the pulsar spin-up rate separately imply relations between these 
two quantities. The size and temperature of the polar-cap hot spot are determined. The 
observed pulse fraction is used to constrain the orientation and geometry of the neutron 
star, and to derive the distance to the system from the calculated luminosity. 



2. Nature of the Stellar Companion 



The present mass M* = M* (i) of the neutron star's stellar companion can be calculated 
as a function of orbital inclination angle i using the measured (Finger, Wilson, & van 
Paradijs 1996b) value of the mass function 
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The solution to equation (1) is plotted in Figure 1, assuming a neutron star mass 
M NS = 1.4M L4 M with M 1A = 1. We find that M*(i) £ 1M unless i £ 6°. The mean 
separation between the neutron star and its companion is given by 



am 



a x sin i 



M NS + M,(i) 
M m (i) sini 



(2) 



where a x sini = 1.12 R Q = 7.76 x 10 10 cm (Finger et al. 1996b). As can be seen in Figure 
1, the values of a(i) range from ~ 25 R Q for i > 5° to ~ 80 R Q for % = 1°. 

Mass transfer in a low-mass system occurs through Roche-lobe overflow rather than via 
a stellar wind (e.g., Verbunt 1990). Given M*(i) and a(i), the average radius of the Roche 
lobe around the stellar companion, i? RL (i), can be determined from the expressions 



R RL (i) ~ a (i) {0.38 + 0.2 log [M,(i)/M NS ]} ; 
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(e. g. Verbunt 1990). -Rrl(0 is plotted in Figure 1 and varies from 4 R Q at i = 90° to 
40 R Q at i = 1°. Such sizes greatly exceed the radii of main sequence stars with masses 
M*(i), which we also plot in Figure 1 (e.g., Mihalas & Binney 1982). Accretion via Roche 
lobe overflow can occur however if the companion star has evolved off the main sequence, 
which is possible within the lifetime of the Galaxy for stars with initial masses ^ 1.0 M & . 
A companion with a current mass > 1 M requires that the system be nearly face-on, with 
% ^ 6°, but mass transfer to the neutron star could have decreased the companion's present 
mass even if its initial mass was > 1 M & . 

The evolved companion star must presently have a radius -R* ^ -RrlW- The radii of 
evolved low-mass stars depend mainly on their helium core masses (Webbink, Rappaport, 
& Savonije 1983). We can use the relations of Webbink et al. (1983) to place a lower limit 
on the present mass of the companion by requiring it to be larger than the helium core 
mass. In the inset to Figure 1, we plot the stellar radius of an evolved companion versus 
its helium core mass as well as the Roche lobe radius versus stellar mass. We find that 
the smallest helium core mass which allows for Roche-lobe overflow in GRO J 1744-28 is 
~ 0.22 M Q . Thus M* ^ 0.22 M and therefore i < 18°. The chance probablity for viewing 
the system with this inclination is 1 in ~20. Also shown in Figure 1 is the no-eclipse 
condition assuming that R* = -Rrl- This restricts i ^ 80°, which is less constraining than 
the above considerations. 



3. Accretion Dynamics 

The Alfven radius for spherical accretion, obtained by balancing accretion and magnetic 
pressure, is given by 

R A = 1.9 x 10 s pt% 7 Ml'J (L 38 R 6 )' 2/7 cm, (4) 

where L38 is the pulsar luminosity in units of 10 38 erg s -1 , fi = B^R^ = 10 30 yU30 G-cm 3 is 
the neutron star magnetic moment {B cq is the surface equatorial field strength), and 10 6 i?6 
cm is the neutron star radius. The corotation radius, where the Keplerian and rigid body 
rotation velocities are equal, is given by R co = (G , M ns ) 1 / 3 (P/2tt) 2 / 3 = 1.0 x 10 s M^ 3 cm 
using P = 0.467 s for GRO J1744-28. The transition from Keplerian motion to corotation 
with the neutron star magnetosphere occurs at r = i]Ra, where 77 « 1. The requirement 
that tjRa < Rco for acretion to take place yields -B cq 3.3 x 10 11 L 1 ^ M{2 Rq 5 ^ 2 r)~ 7 ^ 4 G. 
Requiring magnetic confinement of the hot plasma in the accretion column sets a lower 
limit on the surface polar magnetic field strength B v (= 2B eq for a dipole magnetic field) 
through the expression Bp/8n ^ L/A pc c ~ a^T A /c, where we use the blackbody radiation 
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formula to estimate the polar cap area A pc . Thus B p ^ 2.9 x 10 9 T 2 G, where T 10 is the 
effective temperature in units of 10 keV. 

A more precise estimate for magnetic field strength and accretion luminosity can be 
obtained by interpreting the measured QPO frequency using the beat frequency model 
of Alpar & Shaham (1985), and the pulsar spin-up rate from the Ghosh & Lamb (1979) 
accretion-torque model. Zhang et al. (1996) report a QPO at frequency t'Qpo = 40 Hz with 
integrated power (uncorrected for deadtime) under the 27 Hz FWHM Lorentzian equal to 
5.6% of the total power. Although the width of this QPO is somewhat larger than those 
typically measured in low mass X-ray binaries, with FWHM widths <^ 0.5 the centroid 
frequencies (see, e.g., van der Klis 1995), we still assume that z/q PO can be interpreted as the 
beat frequency between the Keplerian frequency u K = (GM N g) 1 / 2 (27r)~ 1 r~ 3 / 2 at r = r]R\ 
and the neutron star's rotation frequency z/ns (Alpar & Shaham 1985). Less significant 
peaks in the power spectrum are also found at 20 and 60 Hz, and we cannot exclude that 
the 40 Hz peak is a harmonic of the 20 Hz peak but with larger amplitude. This implies 

B eq = 1.1 x 10 10 z>- 7 / 6 L l J 2 M\> 3 Rf' 2 V~ 7/i G (5) 

and 

r]R A = 1.4 x 10 7 i>~ 2/3 M\f 3 cm, (6) 

where we define v = (z^qpo + ^ns)/40 Hz. 

An expression for the period derivative of a disk-fed X-ray pulsar as a function of 
neutron star properties was derived by Ghosh & Lamb (1979). Substituting equation (5) 
into their equation (15) gives 

-P = 2.6 x 10- 5 i)- 1 ' 3 L 38 Mil' 3 Re ^ V~ 1/2 s yr"\ (7) 

where I45 is the neutron star moment of inertia in units of 10 45 g-cm 2 , and f(u s ), proportional 
to the dimensionless torque function (see Ghosh & Lamb 1979; Wang 1995), is a function 
of the fastness parameter u s = (rjR A / R co ) 3 ' 2 . For GRO J1744-28, u s = 0.05 z^ 1 < 1, 
implying f(uj s ) tv 0.9. The observed period derivative is P = —5.95 x 10 -5 s yr _1 (Finger 
et al. 1996b), which could however be different than the value at the epoch 18-19 January 
1996 when RXTE measured the 40 Hz QPO. In terms of the Eddington luminosity 
-^Edd = 1-76 x 10 38 M 1A erg s" 1 , we obtain a pulsar luminosity given by 

L/L Edd = 1.4 77 1/2 *> 1/3 M[ 2 4 /3 Rs 1 J45 • (8) 

Substituting equation (8) into (5), we obtain 

B cq = 1.8 x 10 10 ?T 3/2 M\[l R 6 3 l\l 2 G • (9) 



- 5 - 



The derived value for the magnetic field at the neutron star's equator depends only upon 
quantities intrinsic to the neutron star and details of the accretion process, and is bracketed 
by our estimates given above. We suggest a search with RXTE or ASCA for cyclotron 
absorption features with a fundamental energy at ~ 1-2 keV; note that a larger field is 
inferred if one assumes that the 20 Hz peak in the power spectrum is the fundamental. The 
magnetic field given by equation (9) does not conflict with OSSE observations, which show 
no evidence for cyclotron features at > 40 keV (Strickman et al. 1996). 



4. X-Ray Pulsar Geometry 

Equation (6) for the stopping radius yields a size of the polar-cap hot spot if we 
assume that its area is determined by the footprints of those magnetic field lines that 
cross the magnetic equator at r = i]R A . For an assumed dipolar magnetic field geometry, 
sin 2 m /r = constant, where 6 m is the magnetic colatitude. Thus for an aligned rotator 
9 pc = arcsin 



0.26 z^ 3 M[} /6 R l J 2 



= 15°, although the hot spot size of an oblique 
rotator could be smaller if the matter does not accrete uniformly onto the polar cap. The 
hot-spot size and source luminosity from equation (8) imply an effective temperature of 
= 4.9 r] 1 ^ 8 ^ 1 ^ 12 M^ 2 Rq 3 ^ 4 I^ 4 keV. This is in reasonable agreement with the values of 6 and 
14 keV for the blackbody temperature and e-folding energy reported by BATSE (Briggs et 
al. 1996) and RXTE (Swank et al. 1996), respectively. 

BATSE reports a pulse fraction in the persistent emission of 40% in the 25-45 keV 
band (Finger et al. 1996a), whereas smaller energy-dependent values are measured with 
RXTE for the band between 3 and 12 keV (Zhang 1996). These pulsations are probably 
due to projection effects of the X-ray hot spots on the rotating neutron star. Figure 2 shows 
the various orientations and obliquenesses that produce pulse fractions of 20, 30, and 40% 
calculated from a purely geometric model consisting of two identical uniformly-emitting 
polar caps. The large luminosities radiated from a small area may produce a fan-beam 
radiation pattern, but are not considered here. We include shadowing effects from the 
neutron star but not from the accretion column. A large range of angles are permitted 
for a given pulse fraction. During the bursting episodes, the flux increases by a factor of 
2-8 whereas the spectrum remains nearly identical (a relative excess below 3-4 keV is seen 
during outbursts by RXTE; Zhang 1996). Interpreting the increase in flux as due to an 
increase in emitting area, we predict that the pulse fraction should decrease by an amount 
dependent on a (see inset to Figure 2). The energy-dependence of the pulse fraction could 
be a consequence of thermal stratification, with the hotter inner regions therefore showing 
a larger pulse fraction. 
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The 2-100 keV energy flux reported by RXTE at epoch 18-19 January 1996 (Swank et 
al. 1996) is = 2 x 10~ 7 ergs cmr 2 s _1 for the persistent emission. The luminosity from 
equation (8) then implies a distance to GRO J1744-28 of 

d = 4.6 r// 4 u 1 ' 6 Ml 1 ! 4 /2 b 1/2 R % 1/2 kpc, (10) 

where the geometric factor b depends on neutron-star parameters. By approximating the 
X-ray hot spot as a flat disk, we find b = cos# b s cosa: in the limit # bs + @ pc + ot <C 7r/2 
where neutron shadowing effects are unimportant. For 2-3 magnitudes of extinction per kpc 
in the direction of the Galactic Center, we estimate that the brightest optical counterpart 
would be at V ^ 20 - 24 for an MO III giant {M v = -0.4). 



5. Discussion 

GRO J1744-28 shows characteristics of both X-ray bursters and X-ray pulsars. Its 
intermediate magnetic field strength between the weak fields (B 10 8-9 G) thought to 
reside in low-mass X-ray binaries and the strong fields (B ps 10 12 ~ 13 G) found in X-ray 
pulsars suggests that its burst properties are related to those of X-ray bursters. But the 
relative fluences in the burst and persistent emissions probably rule out a thermonuclear 
origin in the initial bursting phase (Kouveliotou et al. 1996a; Strickman et al. 1996), 
although the energetics derived from RXTE observations (Swank et al. 1996) are favorable 
to such an interpretation. An accretion-disk instability seems implied by the similarity 
of the spectral shapes of the burst and persistent emission, since approximately equal 
energy per particle will be dissipated for an accretion-powered burst. We also note the 
similarities of the burst precursor and post-burst recovery period between the bursts from 
GRO J1744-28 and the Type II bursts of the rapid burster (compare Dotani et al. 1990; 
Lubin et al. 1992), the latter of which are thought to be due to accretion instabilities. 

The beat-frequency model implies a testable relation between the persistent luminosity 
and the QPO centroid frequency given by equation (8). The inferred persistent luminosity 
implies a super-Eddington flux for the bursting episodes, since the burst flux is ~ 2 — 8 
times as great as the persistent emission (Swank et al. 1996). The 3-5 minute recovery 
period following a burst observed by RXTE could be due to radiation effects on the disk. 

The onset of detectable X-ray emission from this source cannot be due to orbital effects 
as in the case of Be X-ray binaries (the system is nearly circular with eccentricity e < 0.026; 
Finger et al. 1996b). We speculate that as the evolved companion expands and overflows 
its Roche lobe, the increased X-ray illumination drives an enhanced fueling episode (see 
Tavani & London 1993) which ceases when the heated atmosphere no longer overflows the 



-7- 

Roche surface. The star then contracts to a quiescent phase after which it resumes its slow 
expansion leading to the onset of another episode of mass transfer and luminous X-ray 
emission. Further calculations of evolution, accretion, and radiation will be required to 
yield a more complete understanding of this extraordinary system. 
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Fig. 1. — Relations between stellar masses and radii which are consistent with orbital 
parameters of GRO J1744-28. We plot the mean separation a, the mass of the stellar 
companion M#, the main sequence radius -R* !ms of a stellar companion of mass M*, and the 
average Roche lobe radius -Rrx as a function of the orbital inclination angle i, assuming a 
1.4 M neutron star. If the companion fills its Roche lobe, the no-eclipse condition requires 
a cos(i) > -Rrl- The inset shows the Roche lobe radius -Rrl (dashed curve) as a function of 
M* and the radius R^ g (solid curve) of an evolved companion as function of its helium core 
mass M c . 
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Fig. 2. — Values of a and 6> obs giving pulse fractions pf = 0.2, 0.3, and 0.4 for X-ray emitting 
hot spots with half-angle 9 pc = 15° (solid curves) and 30° (dashed curve). The obliqueness a 
is the angle between the rotation and magnetic dipole axes, and 6 f, s is the angle between the 
rotation axis and the direction to the line-of-sight. We define pf = (0 max — 0mm) / (0max+0mm)> 
where m ax(mm) represents the maximum (minimum) fluxes in the pulse profile. Inset shows 
variation of pf for polar cap sizes 9 pc = 30° and 60°, using values of a and # b s which yield 
a 40% pulse fraction when 9 pc = 15°. 



